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Abstract 

Using effective theories for jets and heavy quarks it is possible to prove that the double dif- 
ferential top-antitop invariant mass distribution for the process e + e~ — ► ti in the resonance 
region for cm. energies Q much larger than the top mass can factorized into pcrturbatively 
computable hard coefficients and jet functions and a non-perturbative soft function. For in- 
variant mass prescriptions based on hemispheres defined with respect to the thrust axis the 
soft function can be extracted from massless jet event shape distributions. This approach 
allows in principle for top mass determinations without uncertainties from hadronization 
using the reconstruction method and to quantify the top mass scheme dependence of the 
measured top quark mass value. 
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Using effective theories for jets and heavy quarks it is possible to prove that the double 
differential top-antitop invariant mass distribution for the process e + e~ — > ti in the 
resonance region for cm. energies Q much larger than the top mass can factorized into 
perturbatively computable hard coefficients and jet functions and a non-perturbative 
soft function. For invariant mass prescriptions based on hemispheres defined with re- 
spect to the thrust axis the soft function can be extracted from massless jet event shape 
distributions. This approach allows in principle for top mass determinations without 
uncertainties from hadronization using the reconstruction method and to quantify the 
top mass scheme dependence of the measured top quark mass value. 

1 Introduction 

Precise measurements of the top quark mass are among the most important (standard) tasks 
of the ILC project as the top quark mass affects a number of interesting observablcs cither 
directly or indirectly through quantum effects. To be useful such top mass measurements 
have to have small uncertainties, but also need to provide information to which mass scheme 
the measured number refers to. Both aims can be achieved from a threshold scan of the 
cross section cr(e + e~ — * tt) for y/s m 2n%t, from which one expects measurements of the 
threshold masses, such as the IS mass, with uncertainties of about 100 McV [2j j3j 01 [5j [6] . 
Another method is based on mass reconstruction which relies on the idea that the peak 
of the invariant mass distribution of the top decay products is related to the top quark 
mass. This method can be applied at any cm. energy and might also yield uncertainties 
well below 1 GeV [7]. However, until recently it was unknown for which mass scheme 
such measurements can be carried out with small theoretical uncertainties. This is because 
the naive relation between the observable peak of the invariant mass distribution and the 
perturbative top quark propagator pole is affected by hard (i.e. computable) as well as soft 
(i.e. non-perturbative) QCD effects and the present MC tools do not contain the required 
information in a systematic form. Obviously the top mass measurements at the LHC [9] 
suffer from the same problem, but the associated theoretical systematic uncertainty might 
be considerably larger than at the ILC. 

2 Factorization Theorem 

In Ref. |8| a theoretical formalism was presented which remedies this situation, as a first step, 
for the Linear Collider framework, where one does not need to account for QCD radiation 
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Figure 1: Six jet event initiated by a top quark pair, tt — > bWbW — » bqq'bqq'. The plane 
separating the two hemispheres is perpendicular to the thrust axis and intersects the thrust 
axis at the interaction point. The total invariant mass inside each hemisphere is measured. 
Our analysis applies equally well to the lepton+jets and the dilepton channels (not shown). 



arising from the initial state. Assuming a cm. energy Q m t , m t being the top quark 
mass, one can employ the hierarchy of scales 

Q > m > T t > Aq C d (1) 

to establish a factorization theorem for the doubly differential top-antitop invariant mass 
distribution in the peak region around the top resonance: 

M 2 t - - to 2 ~ to T < to 2 . (2) 



dMi dMf 



The invariant masses Mf = (J2 iGXt Pi) 2 ^ = (J2iex { Pi) 2 depend on a prescription X t j 
which associates final state momenta p% to top and antitop invariant masses, respectively. 
For invariant masses in the resonance region the events are characterized by energy deposits 
predominantly contained in two back-to-back regions with opening angles mt /Q associated 
with the energetic jets or leptons from the top decay plus collinear radiation, and by addi- 
tional soft radiation populating the regions between the jets, see Fig. [T] We assume that 
the prescriptions X t i assign all soft radiation to either M 2 or M 2 where the probability 
of radiation being assigned to X t or X% increases to unity when it approaches the top or 
antitop direction. The result for the double differential cross-section in the peak region at 
all orders in a s and to leading order in the power expansion in m t a s /Q, m 2 /Q 2 , T t /m t and 
M t i — m t is given by [5] 



da ( Q 



dMf dM? 



M? - m 2 T 

st.t - — 

mj 



Co H Q (Q,fi m )H m (mj, , Mm, 

V TO,/ / 

dl + dl-B+(s t -—,T ul x) B-(st-—,T t ,Li)s{e + ,r,fi). (3) 

V TO J / V TO J J 

In Eq. ^ the normalization factor cto is the total Born-level cross-section, the Hq and H m 
are perturbative coefficients describing hard effects at the scales Q and toj, B± are pertur- 
bative jet functions that describe the evolution and decay of the the top and antitop close 
to the mass shell, and S is a nonperturbative soft function describing the soft radiation be- 
tween the jets. The result was derived using the hierarchy of scales (p]), matching QCD onto 
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Figure 2: Sequence of effective field theories used to compute the invariant mass distribution. 



(Soft-Collinear Effective Theory) SCET [10] at the scale \x = Q, which in turn is matched 
onto (Heavy Quark Effective Theory) HQET [IT] at a scale /i m of order m t generalized for 
unstable particle effects associated to the large top width T t [12]. An illustration of this 
scheme is shown in Fig. [2j For details on the (admittedly non-trivial) derivation and on 
technical aspects of the factorization theorem we refer to Ref. [5] . In the following we will 
discuss the important ingredients of the factorization theorem and their physical interpreta- 
tion and show what we can learn from them concerning the measurements of the top quark 
mass from the reconstruction method. 



3 Jet Functions and Short-Distance Top Mass 

The coefficients Hq and H m in Eq. ([3]) arise from matching and running in SCET and 
HQET down to the low energy scale /i where one evaluates the jet functions B± and the 
soft function S. These hard coefficients only affect the overall normalization of the invariant 
mass distribution and we will therefore not talk about them here. So let us concentrate 
on the jet and the soft functions, which determine the shape of the distribution and the 
location of the resonance peak. The jet functions describe the perturbative contributions 
of the shape of the invariant mass distribution and are defined by the imaginary part of a 
T-product vacuum matrix element. For the top quark it is 



B + (s t ,T t , n) = Im 



47riV c m,/ 



d A xe ir - x (0|T{^ + (0)W n (0)T^(x)^ + (a;)}|0) 



, (4) 



where v+ is the top four velocity (v+ = 1) and s t = 2v + .r and h v+ is the (HQET) heavy 
top quark field. The vacuum matrix element also contains Wilson lines of the form 

W^(x) = P cxp (ig J dsn - A + (hs + x) S j , W n (x) = P exp ^ — ig J dsn - A + (ns+x) 

where n is a light-like four vector pointing in the antitop direction and A + is field describing 
a gluon that is collinear to the quark. Up to the Wilson lines the vacuum matrix element is 
in fact a heavy quark propagator and, indeed, at tree-level it is just 



Bf e (s,r t ) =Im 



-1 1 
nmj s — 25m + iT t 



1 



i J (s-26m) 2 + r 
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which is the imaginary part of the heavy quark propagator supplemented by a constant 
width term and describing a Breit-Wigner distribution having a width T t . The residual 
mass term 5m becomes relevant at higher orders and controls the mass scheme that is used. 
For the pole mass scheme 5m = to all order in a s . It is the width term (which we can 
approximate as a constant since we are interested in the resonance region) that allows us to 
use perturbation theory for computing the jet function. To understand the role of the Wilson 
lines recall that the two-point function of simple heavy quark fields, evaluated off-shell, is not 
gauge invariant, a fact that becomes e.g. apparent from the gauge parameter dependence 
of the perturbative corrections. The jet functions B±, however, are gauge-invariant due 
to the Wilson lines and well-defined physical objects. Physically the Wilson lines describe 
gluons radiated from the antitop (moving along the four vector n) that are collinear to the 
top quark (moving along i>_|_), and it this additional radiation that renders the jet function 
gauge-invariant and physically meaningful. In momentum space the Wilson lines lead to 
additional Feynman diagrams having l/{n ■ k ± iO) eikonal propagators. 

Having defined the jet function it is now straightforward to address the question which 
mass scheme one might employ to have a good perturbative behavior of the jet function. At 
one-loop [5] one finds that the peak position is located at s pGa k = 25m— CFCt s (fi)/2T t [ln(^-)+ 
|]. Recalling also that the pole mass contains a nasty O(Aqcd) renormalon, it therefore 
natural to use a mass scheme different from the pole mass that is renormalon free and ab- 
sorbs at least the major part of the higher order corrections to the peak position such that 
the resulting scries is convergent. The definition of such a scheme is obviously not unique 
and can also be defined from moments of the distribution [T3|. Gencrically we call such a 
mass a "jet mass" mj and its perturbative relation to the pole mass reads 

??!,/ = Wj 06 — 5m , (5) 

where the HQET power counting requires that 5m ~ a s T in the resonance region. Using 
this relation one can relate the jet mass to other mass schemes. The jet mass mj has already 
been used in the formulae shown before. From this examination we see that top mass one 
can measure from reconstruction is a jet mass. For sure, one cannot measure the MS mass 
from reconstruction because it has 5m ~ ot s mt 3> Tt and would invalidate the HQET power 
counting. 

4 Soft Function 

The soft function S(£ + ,£~,(J,) describes the non-perturbative contributions of the invariant 
mass distribution in the resonance region. Its definition depends on the details of the 
prescription how soft radiation is associated to M t and Mj. One possible prescription is 
using a hemisphere mass definition, where X t and Xi contain everything to the left or 
right of the plane perpendicular to the thrust axis of each event, see Fig. [TJ It is easy 
to understand that such a (and any other) prescription is leading to a non-perturbative 
soft function since one cannot compute perturbatively how the soft particles are distributed 
around the hemisphere boundary. Other prescriptions are possible as long as they do not 
associate soft radiation going in the top direction to the antitop and vice-versa. This is 
in contrast to the jet functions which, according to the condition on X t and Xi stated in 
Sec. [2] are prescription-independent since they describe energetic jets within a small cone 
with opening angle m t /Q around the top direction. At leading order in the power counting 
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the allowed prescriptions do not affect these energetic jets. For the hemisphere prescription 
the soft function is defined by the vacuum matrix element 



s hcmi (e + ,n = w Y<W + - k t a )S(r-K b )(o\(Y n y d (Y n r(o)\x s )(x s \(Y^ (yi) d/ (o)|o) , 



where c, d, e, / are color indices and the Y's are Wilson lines with soft gluons of the form 



Y n (x) = P exp I — ig / dsn ■ A s (ns+x)) , Y^ (x) = P exp I ig / dsn ■ A s (ns+x) 



Yn (%) = P exp ( ig / dsn ■ A s (ns+x) j , Y n (x) = P exp I — ig / dsn • A s (ns+x) J .(6) 



The fc+ a and fcj 6 are operators that pick, according to the hemisphere mass prescription, 
the total + and — light-cone momentum of the gluons that are in hemisphere a and 6, 
respectively, see Fig. [TJ These Wilson lines describe soft radiation off the top and antitop 
quark and also render the soft function gauge-invariant. 

The factorization theorem (J3j> shows that the soft function needs to be convoluted with 
the jet functions. This can be understood physically, since the way how the soft radiation 
is associated to M t and Mf has to affect the observed invariant mass distribution. Field 
theoretically this convolution arises from the fact that the small components of light-cone 
momenta in the top and antitop jets fluctuate at the same length scales as the soft mo- 
menta described by the soft function. At this point it is also useful to note that S is a 
rcnormalized object and that its renormalization group evolution can be computed in per- 
turbation theory. Nevertheless the actual form of the soft function (i.e. the initial condition 
for the soft function evolution at a low energy scale) is not computable with perturbative 
methods. So in practice the soft function needs to be modeled and eventually fixed by ex- 
perimental data, similar to parton-distribution functions. How a soft model function can be 
constructed incorporating consistently the required higher order perturbative information 
has been discussed in Rcf. [14] . 

Given that the soft function is nonperturbative and affects the invariant mass distribution 
at leading order, one might ask what one has gained from predicting the invariant mass 
distribution based on ([3|) and concerning a precise measurement of the top mass from the 
mass M t j where the resonance is located. The crucial aspect is that the soft function is 
universal and appears also in the factorization theorem for event shape distributions for jets 
originating from massless quarks [15j in the dijet region, where the thrust T ss 1 [16] . This 
is related to the fact that the soft Y Wilson lines that arise from massless and from massive 
quark lines are identical. So our factorization theorem for the top invariant mass distribution 
in the resonance region becomes predictive after having determined a soft function from 
event shape distributions from e + e~ data, which arc already available from LEP [17j . Such 
a determination of the soft function was carried out by Korchcmsky and Tafat in Ref. [18j . 

5 Numerical Analysis at LO 

Using the factorization theorem it is straightforward to carry out a simple LO analysis using 
the tree-level result for the jet functions (i.e. one can set dm = 0) and the soft model 
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Figure 3: (a) Plot of the double differential hemisphere invariant mass cross-section 
d 2 a/dM t dMt in units of ia /T 2 for mj = 172, Q = 4.33mj and T t = 1.43 GeV. (b) 
Dependence of the single differential invariant mass distribution as described in the text on 
the cm. energy Q with the same normalization. 



function determined by Korchcmsky and Tafat: 

where Af is a normalization factor. From fits to e + e~ heavy jet mass and thrust LEP data 
they obtained 

a = 2, b =-0.4, A = 0.55 GeV, (8) 

which wc adopt in the following. The analysis illustrates a number of important features 
related to how the predictions by the factorization theorem depend on the cm. energy Q. 

In Fig. (3^ the double differential invariant mass distribution is displayed for the input 
values m,j = 172, Q = 4.33mj and T t = 1.43 GeV. The conspicuous feature of the predicted 
distribution is that the observable resonance peak it shifted toward a higher value by about 
1.5 GeV. This feature is one of the important properties of a invariant mass prescription 
that assigns all soft radiation to the masses. In the factorization theorem it arises from the 
Q/m factor involved in the convolution over the variables £ ± . Intuitively it can be easily 
understood from the fact that the total invariant mass of a fast moving particle with mass 
m plus a soft momentum increases linearly with the soft momentum and the boost factor 
of the massive particle. This feature is also visible in Fig. [3Jd where the single differential 
invariant mass distribution 

da 2 /" A/ "pp- d 2 a 

= - / dM- t , (9) 

dM t r7 Mlowor dM t dM{ K 1 

is plotted. Here the integration interval [Mi owor , M upper ] is twice the size of the measured 
peak mean half width and centered at the peak mass. For the single differential distribution 
one can relate the peak location approximately to the first moment of the soft function by 

Mr^rnj + ^-S^- (10) 

Interestingly this relation can also be used for fixed Q/m for invariant mass prescriptions 
that differ in the treatment of the soft radiation and lead to different first moments of the 
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soft function. When extrapolated to zero moment linearly one can obtain an estimate for 
the jet mass. 

We also find that the distribution gets wider with Q/m. This is again a consequence of 
the Q/m factor occurring in the convolution over £ in the factorization theorem since for 
increasing Q the jet function gets effectively smeared over a wider distribution. The shift 
of the peak position as well as the widening of the invariant mass distribution have been 
observed in simulation studies at the ILC [7] and the LHC for large pr events [§] and can 
now be better quantified using the factorization theorem. 
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